February 2, 2008 9:6 Proceedings Trim Size: 9in x 6in 



AGUIKorea 



ENERGY RELEASE DUE TO ANTINEUTRINO 
UNTRAPPING FROM HOT QUARK STARS * 



D. N. AGUILERA 



t 



Fachbereich Physik, Universitat Rostock, Universitatsplatz 1, 18051 Rostock, 

Germany 

Instituto de Fisica Rosario, Bv. 27 de febrero 210 bis, 2000 Rosario, Argentina 
E-mail: deborahQdarss. mpg. uni-rostock. de 

D. BLASCHKE 

Fachbereich Physik, Universitat Rostock, Universitatsplatz 1, 18051 Rostock, 

Germany 

Bogoliubov Laboratory for Theoretical Physics, JINR, 141980 Dubna, Russia 
email: david@thsunl.jinr.ru 

H. GRIGORIAN 

Fachbereich Physik, Universitat Rostock, Universitatsplatz 1, 18051 Rostock, 

Germany 

Department of Physics, Yerevan State University, Alex Manoogian Str. 1, 
315025 Yerevan, Armenia 
email: hovik@darss. mpg. uni-rostock. de 



An equation of state for 2-flavor quark matter (QM) with diquark condensation un- 
der the conditions for compact stars -/3-equilibrium, charge and color neutrality- is 
presented. Trapped antineutrinos prevent the formation of the diquark condensate 
at moderate densities above a critical value of the antineutrino chemical potential 
fi% . The following consequences are presented: 1) The star develops a 2-phase 
structure (fin,, > u% ): a color superconducting QM core and a normal QM shell. 
2)During the cooling, when the temperature is small enough (T < 1 MeV) the 
antineutrino mean free path becomes larger than the thickness of the normal QM 
shell and the antineutrinos get untrapped in a sudden burst. The energy release is 
estimated as ~ 10 52 erg and an antineutrino pulse is expected to be observed. 
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1. Introduction 

It has been proposed that cold dense quark matter should be in a su- 
perconducting state with the formation of a diquark condensate 1 ' 2 . The 
consequences of the diquark condensation for the configuration and the 
cooling behaviour of compact stars have been broadly studied 3 > 4 ' 5 > 6 and 
the question if this phase can be detected by the signatures still remains 7 . 

Also the engine for the most energetic phenomena in the universe like 
supernova explosions and gamma ray burst does not have a satisfactory 
explanation yet 8 and it has been proposed that the energy involved could 
be related with the occurence of the color supcrdonductivity phase 9,1 °. 

Since the pairing energy gap in quark matter is of the order of the 
Fermi energy, the diquark condensation gives a considerable contribution 
to the equation of state (EoS) that is estimated of the order of (A/n) 2 . 
Disregarding relativistic effects, the total binding energy release in the core 
of a cooling protoneutron star has been estimated as (A/fi) 2 M COIC ~ I0 52 
erg. But, if relativistic effects are considered, the gravitational mass defect 
of the cooling star decreases when diquark condensation is included and 
there is no explosive process 6 possible since the color superconductivity 
transition is second order. 

In this work a new mechanism of releasing the energy in an explosive 
way is presented (for the original idea see 11 ). During the collapse of a pro- 
toneutron star antineutrinos are produced by the /3-processes and remain 
trapped due to the small mean free path. This increases the asymmetry in 
the system and therefore the diquark condensation is inhibited at moderate 
densities. So, a two-phase structure devclopcs in the star: a superconduct- 
ing interior and a sorrounding shell of normal quark matter, the latter 
being opaque to antineutrinos for T > 1 MeV 12 . In the cooling process the 
antineutrino mean free path increases above the size of this normal mat- 
ter shell and an outburst of neutrinos occurs releasing an energy of about 
I0 51 — I0 52 erg. This first order phase transition leads to an explosive 
phenomenon in which a pulse of antineutrinos could be observed. 

1.1. Equation of state for 2- flavour quark matter 

A nonlocal chiral quark model for 2-flavour {u, d} and three color {r, b, g} 
superconducting (2SC) quark matter in the mean field approximation is 
used, for details see 6 ' 13 . The order parameters are the mass gap 4>f and 
the diquark gap A for the chiral and superconducting phase transitions 
respectively. As in 13 , the following chemical potentials are introduced: 
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Mg = {p-u+^d)/2 for quark number, fii = (yu u — Hd)/2 for isospin asymmetry 
and /zg for color charge asymmetry. The deviation in the color space is 
considered /xg <C fi q , so the effect of considering its is neglected. 
The quark thermodynamic potential is expresed as 14 

c\? A 2 

Sl q (<j>, A; n q , /zj, T) + fi„ oc = — + — 



2 

"27 2 



/>GO 

■ / dqq 2 (N c - 2){2E 4> + uj[E^-fi q - m, T] 
Jo 

+uj [E^ - Hg + m, T]+uj [Ej, + Mg - Mi> T]+lu [E^ + n q + T]} 

" 2^ J dqq { E + + E - +UJ ^ ~ ^ 

+ fJ-i , T] + u;[E+ - Hi,T\ + u[E+ + M/ ,T]} (1) 

with 

w[£,T] =TIn[l + exp(-£7/T)] . (2) 

The dispersion relations for the quarks of unpaired and paired colors 
are respectively, 

E^ 2 = q 2 + {m + F 2 (q)<t > f (3) 
Ef = (^ ±A t) 2 + F 4 (g)A 2 (4) 

The interaction between the quarks is implemented via a Gaussian form- 
factor function F(q) in the momentum space (Gaussian types give stable 
hybrid configurations 7 ) as F(q) = cxp(— q 2 / A 2 ) . 

The parameters A = 1.025 GeV, G\ = 3.761 A 2 and m u = rrid = m — 
2.41 MeV are fixed by the pion mass, pion decay constant and the con- 
stituent quark mass at T = zx = 15 . The constant G2 is a free parameter 
of the approach and is fixed as G 2 = 0.86 G\. 



1.1.1. Stellar matter conditions 

The stellar matter in the quark core of compact stars is considered to 
consists of u and d quarks and leptons (electrons e~ and antineutrinos 
9 e ) under the following conditions 

• /3-equilibrium d < — ► u + e~ + P e , fi e + fi^ e = — 2/zj, 

• Charge neutrality |n„ — \nd — n e = 0, ng + n/ - 2n e = 0, 

• Color neutrality n s = 0, 2n qr — n q b = 0, 
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where the number densities corresponding to the 

3 T,0=0o,A=A o 

chemical potential Hj defined above. 

The lepton contributions (I = e, v e ) as ideal Fermi gases 

n , ( ,,r) = -A-y-l^-^T- (5 ) 

are added to the quark thermodynamical potential 
0(0, A; ft q , iij, Me , T) = O 9 (0, A; n q , M/ , T) + O e (,i e , T) + Pe (/^ , T). (6) 

The baryon chemical potential [Ib = — ^/ is introduced as the 
conjugate of the baryon number density ns- 

The n function can have several minima in the 4>, A plane, an example 
is shown is Fig. 1. The global minimum represents the stable equilibrium 
of the system and the minima search is perfomed solving the gap equations 



on 



dK 



= 



(7) 



p=<p :A=A UL * 0=0 o ;A=A o 

under the conditions that are mentioned above for the stellar interior. 
The thermodynamics of the system, e.g. pressure P, energy density e, 
number density n and entropy density s, is defined via this global minimum 



0(</>o, A ; /J. B ,Vi,T) = e-Ts- [i B n B - Vinj = -P . (8) 

To fulfill the charge neutrality condition (see Fig. 2, right) a mixed 
phase between a subphase with diquark condensation (subscript A > 0) 
and normal quark matter subphase (subscript A = 0) is defined via the 
Glcndenning construction. The Gibbs condition for equilibrium at fixed T 
and hb is that the pressure of the subphases should be the same 



P = P A>0 (HB,Vl^e,T)=P A =°(LlB,Hl^e,T) . (9) 

The volume fraction \ that is occupied by the subphase with diquark 
condensation is defined by the charge Q in the subphases, 

X = Qa>o/ (<3a>o - <2a=o) , (10) 

and is plotted on the right panel in Fig. 2 for different antincutrino 
chemical potentials as a function of /i^. In the same way, the number 
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Figure 1. Cuts of the thermodynamic potential Q((f>, A; /ig , fii, T = 0) in the planes 
A = const (on the left) and <j> = const (on the right) for two different constant values 
of fig and the corresponding fi[. For fig = 933 MeV (upper panel) two degenerate 
minima can coexist at the values: <f> = 331 MeV, A = (solid lines) and 4> = 107 MeV, 
A = 98 MeV (dashed lines). For fig = 1100 MeV (lower panel) the minimum with a 
nonvanishing diquark A = 121 MeV and cf> = 54.8 MeV (dashed lines) is preferable. This 
corresponds to a first order transition from the vacuum to a superconducting phase. In 
this example G2/G1 = 1 was taken. 

densities for the different particle species j and the energy density are given 
by 

e = xe A >o + (1 - x)eA=o • (12) 
1.1.2. Equation of state with trapped antineutrinos 

Increasing the antineutrino chemical potential fi p<! increases the asymme- 
try in the system and this shifts the onset of the superconducting phase 
transition to higher densities. 

Above a critical value of /i Pe > [i% ~ 30 MeV (see critical value in Fig. 
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Figure 2. Left: Solutions of the gap equations and the charge neutrality condition 
(solid black line) in the fij vs. [i e plane. Two branches are shown: states with diquark 
condensation on the upper right (A > 0) and states from normal quark matter (A = 0) 
on the lower left. The plateau in between corresponds to a mixed phase. The lines 
for the /3-equilibrium condition are also shown (solid and dashed red lines) for different 
values of the (anti)-neutrino chemical potential. The stellar matter should satisfy both 
conditions (intersection of the corresponding lines) and therefore for = a mixed 
phase is preferable. Right: Volume fraction \ of the phase with nonvanishing diquark 
condensate obtained by a Glendenning construction of a charge-neutral mixed phase. 
Results are shown for two different values of /j,a e . 

2, on the left) first a normal quark matter phase occurs and then the phase 
transition to superconducting matter takes place, see Fig. 3, on the left. 

The consequences for the equation of state can be seen on the right of 
the Fig. 3: the onset of the superconductivity in quark matter is shifted to 
higher densities and the equation of state becames harder. 



1.2. Quark stars and antineutrino trapping 

The configurations for the quark stars are obtained by solving the Tolman- 
Oppenheimer-Volkoff equations for a set of central quark number densities 
n q for which the stars are stable. 

In Fig. 4 the configurations for different antineutrino chemical poten- 
tials are shown. The equations of state with trapped antineutrinos are 
softer and therefore this allows more compact configurations. The presence 
of antineutrinos tends to increase the mass of the star for a given central 
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H B [MeV] ^ B [MeV] £[M eV/fm 3 ] 

Figure 3. Left: Solutions of the gap equations and /ij as a function of /ib- Increasing 
the antineutrino chemical potential increases the asymmetry in the system and the su- 
perconducting phase is inhibited at moderates densities. Right: Equation of state for 
different values of of trapped antineutrinos. As /ip c increases the equation of state 
becomes harder. 

density. 

A reference configuration with total baryon number Nb — 1.51 -ZVq 
(where jV© is the total baryon number of the sun) is chosen and the case 
with (configurations A and B in Fig. 4) and without antineutrinos (/ 
in Fig. 4) are compared. A mass defect can be calculated between the 
configurations with and without trapped antineutrinos at constant total 
baryon number and the result is shown on the right panel of Fig. 4) . The 
mass defect could be interpreted as an energy release if the configurations 
A, B with antineutrinos are initial states and the configuration / without 
them is the final state of a protoneutron star evolution. 

1.2.1. Protoneutron star evolution with antineutrino trapping 

After the collapse of a protoneutron star the star cools down by surface 
emission of photons and antineutrinos. Antineutrinos are trapped because 
they were generated by the direct /3-process in the hot and dense matter 
and could not escape due to their small mean free path. The region of the 
star where the temperature falls below the density dependent critical value 
for diquark condensation, will transform to the color superconducting state 
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Figure 4. Left: Quark star configurations for different antineutrino chemical potentials 
= 0, 100, 150 MeV. The total mass M in solar masses (M Bun = Mq in the text) is 
shown as a function of the radius R (left panel) and of the central number density n q in 
units of the nuclear saturation density no (right panel). Asterisks denote two different 
sets of configurations (A,B,f) and (A',B',f) with a fixed total baryon number of the 
set. Right: Mass defect AAf and corresponding energy release AE due to antineutrino 
untrapping as a function of the mass of the final state Mf. The shaded region is defined 
by the estimates for the upper and lower limits of the antineutrino chemical potential in 
the initial state fin^ = 150 MeV (dashed-dotted line) and = 100 MeV (dashed line), 
respectively. 



which is almost transparent to (anti)neutrinos. But nevertheless due to 
the trapped antineutrinos there is a dilute normal quark matter shell which 
prevents neutrino escape from the superconducting bulk of the star, see Fig. 
5 and Fig. 6. The criterion for the neutrino untrapping transition is to cool 
the star below a temperature of about 1 MeV when the mean free path of 
neutrinos becomes larger than the shell radius 16 . If at this temperature the 
antineutrino chemical potential is still large then the neutrinos can escape 
in a sudden outburst. If it is small then there will be only a gradual increase 
in the luminosity. An estimate for the possible release of energy within the 
outburst scenario can be given via the mass defect defined in the previous 
subsection between an initial configuration with trapped neutrinos (state 
A or B) and a final configuration without neutrinos (state /). 
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T = 40MeV T<10MeV 




Figure 5. Left graph: Quark star cooling by antineutrino and photon emission from the 
surface. Middle graph: Two-phase structure developes due to the trapped antineutrinos: 
a normal quark matter shell and a superconducting interior. Right graph: Antineutrino 
untrapping and burst-type release of energy. 
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Figure 6. Star evolution corresponding to Fig. 5 plotted in the phase diagram 



2. Conclusions 

The effects of trapped antineutrinos on the diquark condensates in quark 
star configurations are investigated. At fixed baryon number the energy 
release in the antineutrino untrapping transition is of the order of 10 52 erg. 
This is a first order transition and leads to an explosive release of energy 
that could help to explain energetic phenomena in the universe like gamma 
ray bursts or supernova explosions. 
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